A simple model for magnetic inelastic dark matter (MiDM) 

Sudhanwa Patra^'^'* and Soumya Rao^'''^ 

^Physical Research Laboratory, Ahmedahad 380009, India 
^Institute of Physics, Bhubaneswar 751005, India 

Abstract 

A simple model for magnetic inelastic dark matter (MiDM), which is a minimal extension of the 
standard model with right-handed neutrinos (Nji), a singly charged scalar (S) and a vector-like 
charged fermion (E), has been presented. In this framework in which dark matter inelastically 
scatters off nuclei through a magnetic dipole interaction by making a transition to a nearly degen- 
erate state with a mass splitting {6) of 100 keV. The model is constrained from Direct detection 
experiments, DAMA and COGENT. In our analysis we do not find any parameter space using the 
recent annual modulation data from COGENT, where DAMA and COGENT are consistent with 
each other. The bound from relic density provides a much stronger constraint. We find that this 
minimal extension of SM incorporates a DM candidate which can explain the indirect detection 
results while being consistent with relic density measurements. The right handed neutrino and the 
charged particles responsible for the right-handed neutrino magnetic moments could be produced 
at the Large Hadron Collider. 
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I. INTRODUCTION 



The quest for the identification of dark matter (DM), together with the comprehension 
of the nature of dark energy, is one of the most challenging problems in the understanding 
of the physical world. A canonical model for Dark Matter (DM) utilizing the existence of 
Weakly Interacting Massive Particles (WIMPs) has emerged. On the other hand, supersym- 
metric extensions of the Standard Model naturally incorporate such stable particles, like for 
example the neutralino. The observed DM density, set by thermal freeze-out, determines 
the cross-section to annihilate to Standard Model (SM) fields to be a value typical of weak 
scale physics (av) ~ 3.6 x lO^^^cm^/sec. Within the paradigm of these models, many phe- 
nomenological expectations have been fixed, including the annihilation modes to the SM 
interaction channels with corresponding rates for indirect detection in the galaxy today. 

Cosmological observations strongly suggest that non-luminous, non-baryonic matter con- 
stitutes most of the dark matter in the universe. This dark matter should be distributed 
in dark halos of galaxies such as the Milky Way, enabling the direct detection of the dark 
matter particles via their interactions in terrestrial detectors. Recently, inelastic dark mat- 
ter (IDM) is a viable scenario to explain DAMA consistent with the other experiments. 
The inelastic scenario assumes that WIMPs (x) can only scatter off baryonic matter (N) by 
transition into an excited state at a certain energy above the ground state {xN — > x*^)^ 
while elastic scattering is forbidden or highly suppressed. 

Since neutrinos have very weak interactions, any additional interaction, like the one pro- 
vided by magnetic moments, could have dramatic consequences in their behaviour. As a 
result, neutrino magnetic moments provide new photon-neutrino couplings which affect the 
production and detection of neutrinos at colliders and large transition magnetic moments 
could affect neutrino oscillation. The magnetic moments of right-handed neutrinos arise as 
dimension 5-operators which is the same dimension as of the well-known Weinberg opera- 
tor that is often used to parametrize neutrino masses. It is therefore natural to consider 
right-handed neutrino magnetic moments as the first manifestation of non-trivial electro- 
magnetic properties of neutrinos. Among important low-energy properties of WIMPs are 
their electromagnetic form factors. It has been known for a long time that the possibility 
of charged WIMPs is strongly disfavored [1, 2] and stringent limits exist in the case of a 
fractional charge [3]. 
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Motivated by the discussion above we present here a very simple model which gives 
rise to right-handed neutrino magnetic moments; it includes, in addition to the SM fields 
and the right-handed neutrinos, a charged scalar singlet and a charged singlet vector-like 
fermion. While the phenomenological effects of an effective magnetic coupling of right- 
handed neutrinos have their own intrinsic interests, their observation would also indicate 
the presence of physics beyond the SM, and it is therefore natural to determine the types 
of new particles and interactions that would be implied by the observation of such effective 
couplings. There is, of course, a certain amount of freedom in constructing models based on 
the single constraint that Majorana magnetic couplings be generated at the one loop level, 
but the simplest possibility, involving the introduction of one charged heavy scalar and one 
charged heavy fermion, involves the basic features of more general models and addresses 
the main problems such theories face. In addition, if the right-handed neutrino magnetic 
moments are large the new particles should be relatively light and could be produced and 
analysed at the Large Hadron Collider (LHC) or the International Linear Collider (ILC). 
This is also a general feature of models with large right-handed neutrino magnetic moments. 

Recently there has been a lot of excitement in the field of Dark Matter with the obser- 
vation of an annual modulation signal by COGENT [ i], similar to what has been reported 
consistently by DAMA [ ] over several years. An excess in cosmic ray positron and electron 
signals over the expected background as observed by PAMELA [6] and FERMI [ ] may 
be a signal of annihilating Dark matter for the energy of (9(100) GeV. The annihilation 
cross-section needed to produce these signals is non-thermal, a factor of ~ 10 — 1000 (de- 
pending on DM mass and astrophysical boost factor) large than the thermal annihilation 
cross-section [8, 9]. One possibility is to add new particles and these particles mediate a 
Sommefeld enhancement, implying boosted annihilation in the halo today, while also acting 
as intermediate final states, thereby allowing SM particles produced from DM annihilations. 

The paper is organised as follows. In section-II, we present a simple model implement- 
ing inelastic dipolar dark matter which is a simple extension of the Standard model for 
the dipolar dark matter to satisfy both obsevational and experimental bounds. In the 
following section-Ill, we introduce the effective lagrangian for the magnetic inelastic dark 
matter(MiDM) interaction with a photons. We discuss how one can large magnetic dipole 
moment for heavy Majorana neutrino which is assumed to be dark matter candiadate in our 
model and analytical expressions are presented in section III. We discuss the relic abundance 
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in section IV. Section V presents constraints on dark matter diple moments and masses that 
arises from direct search detectors DAMA and COGENT. Comparison with the indirect 
detection experiments PAMELA and FERMI is discussed in section VI. Finally, we present 
our summary of our whole work in section VII. 

II. MODEL 

Beside the standard model particles and right-handed Majorana neutrinos, the model 
contains a singly charged scalar and three vectorlike singly charged fermions. The new 
fermions are assumed to be vector like to make sure that the theory is anomaly free as for self 
consistency. The quantum numbers of these extra particles under SU (3)c x SU {2)l x U{1)y 
gauge group : Nji = (1,1,0), S = (1,1,0) and E = (1,1,0) and their Z2 assignment is 
shown in Table: [I]. The hypercharge assignment is done by the relation Y = Q — I^. 



TABLE I: Particle content of the proposed Model 





Field 


SU{3)c X SU{2)l X U{1)y 


Z2 


Fermions 


Ql = {u,d)l 


(3, 2, 1/6) 


+ 




UR 


(3, 1, 2/3) 


+ 




dR 


(3, 1, -1/3) 


+ 




iL = (^^, e)l 


(1, 2, -1/2) 


+ 




e-R 


(1, 1,-1) 


+ 




E 


(1, 1,-1) 






Nr 


(1, 1, 0) 




Scalar s 




(1, 2, +1/2) 


+ 




s+ 


(1, 1, +1) 


+ 



We first analyze the case when the DM candidate is a fermion we denote by x, which is 
assumed to be odd under Z2 , and a gauge singlet; we will also assume that x has no chiral 
interactions. In this model, the Yukawa terms of the lepton sector would be 
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c = - 



Ys {NrYES^ + Y's NnES+ + h.c. 



-(iV^)^M^iV^ - MeEIEl + h.c. 



(1) 



where Ys is the new Yukawa couphngs, and Me are the mass matrices of iV^ and Ei. 
In the unitary gauge the Higgs doublet is given by 



H + v 10 
^ 1 



where H is the SM Higgs boson and v = 246 GeV is the vacuum expectation value. All 
through this study we will choose Mh = 120 GeV. 

The general scalar potential contains the quadratic and quartic terms as below 



S) = -ml ($t $) + Ai (<l>t {S^ S) + A2 {S^ 

+ ^(5t5)($t$) (2) 

where the mass of the singly charged scalar is m| = + A2f^. 

III. MAGNETIC INELASTIC DARK MATTER (MIDM) 

In this section, we focus on magnetic inelastic dark matter (MiDM), because it has a 
unique and interesting directional signature and it has been shown that MiDM could explain 
both DAMA and other null results [10]. The model takes advantage of both the magnetic 
moment and large mass of iodine. In [I I] the velocity dependent scattering resulting from 
magnetic dipole operator is said to explain simultaneously DAMA and COGENT results 
as well as null results from other experiments. Explicit formula for scattering cross section 
assuming magnetic dipole operator has been derived in [12]. 

In MiDM, the dark matter couples off-diagonally to the photon: 

j^^-^r^.uJ^-'x (3) 
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where mass of x ^i-nd x* ^"^^ split hy 5 = 100 keV. The off-diagonal coupling is natural if the 
dark matter is a Majorana fermion. In the present scenario, heavy Majorana neutrino is a 
candiadate for magnetic inelastic dark matter (MiDM). Due to the Majorana nature, the 
magnetic moment of heavy Majorana neutrinos is zero. There is only transition magnetic 
moment for them. The effective Lagrangian for coupling of a heavy Majorana neutrino with 
a magnetic dipole moment /z 

where /i^-fc is the transition magnetic dipole moment. 




FIG. 1: Contributing diagram for the heavy right handed Majorana neutrino 
A. Magnetic moment of heavy RH Majorana neutrino 

In the model considered, we have four diagrams contributing to the transization magnetic 
moment, which are depicted in Fig.(l): (a and b) a loop with the B gauge boson attached 
to the E and (c and d) a loop with the B gauge boson attached to the scalar S. 

The Yukawa interactions of heavy Majorana neutrinos with x ^-iid can be rewritten in 
the following way 
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2 L 



2 L 



tC 



(5) 



through which we can derive relevant feynman rules. 

Assuming that heavy Majorana neutrinos are nearly degenerate, i.e., Mr, 
we derive the expression of /u(0): 



[X{MImIMe,)-X{Me,,MImI)\ , (6) 
with 

X(A, B,C) = j + + ' 

where Me^ and Ms are the mass eigenvalues of heavy vector-like fermion E and singly 
charged scalar 5, respectively. In simplest form, one can write as 

/^./.(O) ~ ^,{Yl)p,{Ysho. ^F{x) (7) 

where F{x) = + ^j3^1n(a;), x = Non-perturbatibe limit gives us (YgYs) < 47r. 
With this spectrum, one can get the magnetic moment of the order of 10~^/is to lO^^/i^- 



/^a/3(0) 



fi47r2 



('^5)/3il^5; 



ia 



B. Cross-section RH neutrino through e^e — )■ 7, Z* — )^ N^Nj [k / j)] 

The most dramatic effect of a large magnetic dipole moment of a heavy neutrino will 
be in the production cross section and angular distribution. Though the discussion of the 
differential cross section for a heavy charged lepton can be found in work of Sher [13], here 
we only discuss qualitatively how one can produce RH Majorana neutrinos in near future 
experiment within our framework. In the discussion of Escribano and Masso [1 I], one can 
write a U{1) invariant operator as: Nrj {^^^ + iT)-'^) aap N'j^f^B"'!^ , where 5"^ is the U{1) 
field tensor. This gives a coupling to the photon, which we define to be the EDM, as well 
as a coupling to the Z which is the EDM times ianOw When we include the effect of Z 
coupling to in the differential cross section, it turns out that the contribution has very 
little effect on the result. 
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The differential cross-section for the process, e~^e — )■ 7, — )■ NkNj {k 7^ j), is given by 

'1- — ( + , . PzzJ'2 



da 

dn 



a' 
4^ 



+ 



'1—4 sin^ 9w) tan 9w 



sin^ 29 



w 



8 sin^ 2^11/ 



where the values of Fi, J-2, J-3, P^z and P^z 



= ix%s sin^ 6* ( 1 + 



4M^ 



= 1 + cos^ 



4M2 



sin^ 9 + 8Cv cos 9 



+ filf s tan^ 9 



w 



S\Yi^9+ (l + cos^^) 

s 



J^3 = s 



sin^ ^ + (1 + cos^ 9) 

s ^ ' 



P. 



zz 



7Z 



[s - M|)2 + r2M| 
b{b - Ml) 



with /i^ is the transition magnetic moment of heavy Majorana neutrino, C\ 
and we have dropped the numerically negligible Cy terms, for simplicity. 



(8) 



(9) 

-2 sin^ 6'vK, 



C. Mass splitting between two Majorana neutrino component 

In this paper we shall not discuss the origin of the small mass differences between the 
degenerate right-handed neutrinos, but for completeness we demonstrate that a mass split- 
ting of the order of 10^'' GeV is not completely unnatural for 100 GeV scale right-handed 
neutrinos. Consider a diagram with a vertex A3($^$)(S'^S') attached to the singly charged 
scalar S which runs in loop and this kind of diagram gives a finite contribution to the mass 
splitting. The diagram with a vertex A3($^$)(S'^5') attached to the 5* field (which is shown 
in fig: (3) give a finite contribution to the mass splitting as: 

AM«~M1^<?>^ (10) 

^ (47r)2 ^ME ^ ' 

For the mass of the charged lepton to be around 100 GeV (i.e. Me ~ 100-150 GeV), ($)=100 



GeV and Y| ~ 1 — 10, one can write 



10^ 

AMr AccE, (11) 
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FIG. 2: The total cross section for heavy right handed neutrino e^e~ — )• 7, Z* — )• N^Nj (k 7^ j) 
for various EDMs, in units of Bohr magneton. The cross section is shown as a function of center 
of colhder energy y/s and here we have varied the masses of heavy right handed neutrino as 
M = 100, 160, 250, 400, 500 GeV from the top to the bottom curves. 
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FIG. 3: Mass splitting between the dipolar dark matter(DDM) and its companion partner 



Now one can easily get the mass splitting between two right handed neutrinos of the order 
of O(10~^) GeV. If we thus start with a symmetry to get a 100 GeV Scale degenerate 
right-handed neutrinos, after the symmetry breaking, we get a mass splitting between the 
companion states of right-handed neutrinos to be in the range of 0(10^^) GeV or 100 keV 
range, naturally through radiative corrections. 
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IV. DARK MATTER ANNIHILATION AND RELIC ABUNDANCE 



To calculate the relic abundance we solve the Boltzmann equation for heavy Majorana 
neutrino which is assumed to be a Dipolar Dark Matter candidate, given by [1 ">] 

^ + 3Hn^ = -{av) K " ^.cj 
where H is the Hubble parameter, is the number density of Dipolar Dark Matter (x), 
n^gq is the equilibrium number density, v is the relative velocity and (cv) is the thermal 
average of the annihilation cross-section xx ~^ ^ill. 

Since the Dipolar Dark Matter considered in our model is non-relativistic, the equilibrium 
number density is n^,cq = 9xi^xl (27rx))^/^ produced by the back-reaction //' — )• % x at 
a given temperature, where is the spin degrees of freedom, T is the temperature and is 
the mass of the relic. For particles which may potentially play the role of cold dark matter, 
the relevant temperature is order of ^ and the non-relavistic equilibrium abundance is well 
justified. 

In terms of the dimensionless variables = n^/s where s is the total entropy density, 
and X = ^ the boltzmann equation becomes 

rfl^ ^ _A(X) r 2 _ y2 1 /-,2^ 

dx ^ x.cqJ v-^ / 

where 

\(x) ^ (^) Af„,M, (av}(x) (13) 

n«-Sfl)">^''V- (14) 



2tt^ \8J 

where g^, and g^:s are the effective degrees of freedom of the energy density and entropy 
density respectively. Here Mpi = 2.4 x 10^^ GeV is the Planck mass and is the mass of 
Dipolar Dark Matter (x). 

In the case of Majorana fermions only non-identical fermions can annihilate {X1X2 ff) 
through the dipole moment operator as only the transition dipole moments {'D12, fii2) are 
non-zero. When the mass difference 6 between xi aiid X2 is small, the cross section for 
Majorana annihilation process is identical to that of the Dirac fermions (with P, /i replaced 
with T>i2,fii2)- The iDDM pairs annihilate to either photons or charged pairs through the 
diagrams shown in Fig: 
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FIG. 4: Annihilation diagram for dipolar dark matter(DDM) pair into fermion particle-antiparticle 
pair 



Hence the total cross-section is crann"^ = o'xx*->27'^ + '^ff^2y'^- Here we will calculate the 
comsological relic abundance Qh'^ of the iDDM by asuuming standard freeze-out of annihi- 
lations via the dipole coupling to 7 and there is no XX asymmetry. In the nonrelativistic 
limit, the thermal average annihilation cross section {(tv) reduces to an average over a 
Maxwell-Boltzmann distribution function \' ']: 



™3/2 fl 



{(^v) = {(xv)v e 4 dv (15) 

The DM annihilation rate for magnetic dipolar interaction, with fermions as the annihilation 
products is given by [16] 

(-) = # (16) 

However we find that the dominant contribution to relic density comes from the annihi- 
lation through magnitic dipole interaction into W^W~ given by 



3eVi2 L _ ^ 
Atx \ m? 



H = ^(6-^) (17) 



X 



The present day mass density of iDDM particles x is given by 

2 Xf 8.5 X 10-1^ GeV-^ 

where we take Xf = 20. 

The observed value of CDM density from the seven year WMAP data is ^cdm h"^ = 
0.1123±0.0105 (3cr) [ ] (where h is the Hubble parameter in units of 100 km/s/Mpc). The 
annihilation cross-section required to get this value of fi/i^ is 1.81 x lO^^^cm'^s^^. To obtain 
this value of annihilation cross section in our model we choose a DM mass of 160 GeV which 
gives = 3.8 X 10^"^ fiB- 
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V. DIRECT DETECTION OF MIDM 



In direct detection of DM one measures WIMPs scattering off target nuclei. The scat- 
tering can be assumed to be either elastic or inelastic. In the case of inelastic scattering, 
the process is xi + ^ — X2 + N where xi and X2 are two different mass eigenstates, and in 
general, there is a mass difference between xi and X2-, ^ = ""^2 — ^i- Due to this mass dif- 
ference, the minimum DM kinetic energy needed for the nucleon scattering becomes higher 
[18]. There is a minimal velocity required to produce recoil energy in such an inelastic 
scattering 

where m^v is the mass of the target nucleus, /i is the reduced mass of the WIMP-nucleus 
system, and 5 is the WIMP-mass splitting. Here we assume the WIMPs to have a Maxwellian 
velocity distribution as described in [19]. 




FIG. 5: Inelastic scattering of MiDM particle with a nucleus through the magnetic moment of 
fermion singlet induced by a charged lepton loop and a photon vertex at 1-loop 

The differential cross section per unit energy transfer for elastic scattering by a magnetic 
dipole moment interaction is given by [12] 



da 



dER A-kEr 



E, 



E, 



+ 



31 



2mr) 



\F{q)\' 



(20) 



where the first term is due to interaction with the charge of the nucleus and the second term 
with the magnetic dipole moment of the nucleus. Here / denotes the nuclear spin, fi^ and 
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jj^ denote the magnetic moments of nucleus and DM respectively, Z is atomic number of 
the target nucleus and is the nuclear form factor. Although the nuclear charge form 

factor and the nuclear magnetic moment form factor can be different we assume them to be 
approximately equal for simplicity [ ]. The differential cross section given by 20 is used to 
calculate the event rate spectrum. We use the recently released modulation data from the 
COGENT experiment [1] as quoted in []' '] and also the latest modulation data from DAMA 
[5] . The allowed parameter space extracted assuming magnetic dipolar interaction is shown 
in Fig. 6. It can be seen from the figure that DAMA and COGENT can both be explained 
for a small sliver of parameter space for a very small mass difference of about 4-5 keV. Also 
the bound from relic density is a much more stringent constraint compared to the direct 
detection bound from DAMA or COGENT. 




-6.2- 



-6.4^ ^- 

50 100 150 200 



FIG. 6: Allowed parameter space in the plane for COGENT and DAMA. Solid line represents 
bound from relic density. 
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VI. INDIRECT DETECTION OF MIDM 



The recently reported excesses of positrons and antiprotons in cosmic rays by the ex- 
periments hke PAMELA [6, 25, 26] and of electrons+ positrons by Fermi [7, 27, 28] can 
be explained by dark matter annihilation with possibly a boost factor required in the DM 
annihilation cross section. This boost factor is explained either by astrophysical sources or 
effects like Sommerfeld enhancement. 

We fix the DM mass to be 160 GeV which gives the correct relic density for a DM 
annihilation cross section of 1.81 x lO^^^cm'^s"^. The dominant contribution to relic density 
comes from the W channel as shown in the previous section. But in order to obtain the 
positron and antiproton fluxes for the indirect detection analysis, we take contributions 
from W as well as the leptonic channels e, /x, and r. In case of e the cross section will be a 
^-function. 

The method followed here to obtain the positron and antiproton fluxes is outlined in 
[21]. We use [22] to obtain event spectra for positrons and antiprotons produced from iJ 
{I = e, /i, r) and W+W' . This is then fed into GALPROP [23, 21] which calculates the 
final fluxes of positrons and antiprotons to be compared with experiment. The diffusion 
parameters used in the Galprop code are as tabled in [ ' ' ] . The DM annihilation cross 
section used for calculating the positron and antiproton fluxes is 1.81 x lO^^^cm'^s"^, a 
boost factor of 10 is required to fit the data. Our results are shown in Fig. 7, Fig. 8 and 
Fig. 9. The positron flux shows fairly good agreement with the PAMELA data [6] as seen 
in Fig. 8. The antiproton flux is also within the observed values as shown in Fig. 9 due to 
a relatively light DM mass of 160 GeV and a very small boost in the cross section. 

VII. SUMMARY 

We have consider a minimal extension of the standard model which can explain non-zero 
light neutrino mass, magnetic moment of heavy RH singlet neutrino which we consider here 
as a DM candidate. We have shown that the inelastic dipole dark matter particle (which 
is a heavy RH singlet fermion in this case) with non-zero electric and/ or magnetic dipole 
moment can satisfy the experimental and observational bounds. Furthermore, this scenario 
may be tested at future particle colliders (such as the Large Hadron Collider (LHC) ) or 
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FIG. 7: The (e" + e+) flux for the 160 GeV DM compared with FERMI-LAT data [27, 28]. Dashed 
denotes the CR background and dotted hne is the DM annihilation signal. 
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FIG. 8: Positron flux ratio for the 160 GeV DM compared with Pamela data [6]. The annihilation 
cross section is taken to be avrei = 1-81 x lO^^^cm^s"^ which corresponds to a boost of 10. 

dark matter detection experiments. 

At the same time, it gives the connection between neutrino mass and magnetic dipole 
coupling which has opened up new possibilities for model builders. We have considered the 
dipole moment interactions between the heavy-light and heavy-heavy neutrino counterparts. 
As a result, the large magnetic moments of heavy Majorana neutrinos can enhance the 
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FIG. 9: Antiproton/Proton flux ratio for the 160 GeV DM compared with Pamela data [25, 26]. 
The annihilation cross section is taken to be avrei = 1-81 x Id^^^cm^s^^ for all annihilation channels 
with a boost of 10. 

production cross section of TeV scale right-handed neutrinos though the Drell-Yan process, 
e^e~ — )■ 7, — !■ NiNj {i ^ j), which is within the reach of the future linear collider (ILC). 

The model also provides a good agreement with the indirect detection experiments like 
PAMELA and FERMI consistent with the relic density bound from WMAP. At the same 
time we also constrain magnetic dipole moment from Direct detection experiments like 
DAMA and COGENT. In our model we do not find a common parameter space for DAMA 
and COGENT using the annual modulation data from both experiments. We also find that 
the relic density constraint is much more stringent compared to constraints from DAMA 



and COGENT. 
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